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Abstract. It is demonstrated that a certain integral equation can be solved using the 
'] Painleve equation of third kind. Inversely, a special solution of this Painleve equation 

^0 ' can be expressed as the ratio of two infinite series of spheroidal functions with known 

. coefficients. 

Ph- 1. Introduction 

We consider an integral equation in a finite interval for unknown function g(t) 



J K(\x-t\)g(t)dt=f(x) (1) 



where kernel K(w) has the form 



> 
00 

l> 

^. : K(w) = w v K u (6w) . (2) 

Here K v {w) is the modified Bessel function of the third kind which obeys the equation 

K{w) + -K(w) - f i + 4) = ( 3 ) 

c3 I and exponentially decreases at large argument 

K v {w) -> ,/^e- . (4) 



At small w 

KJw) — ► 2^ 1 I»i/;- I/ + 2- v - 1 r(-z/)^ + 0(w 2 - M ) . (5) 

We restrict the discussion to the case \u\ < 1/2 though certain results could be 
generalized to all v > —1/2. 

Equation (pQ) with kernel appeared in some physical problems. For example, it 
describes forced convection heat transfer and contact problems in elasticity (see [TJ [2] 
and references therein). The particular case of ([[]) with v — and imaginary 9 gives the 
solution of the wave scattering problem with the Dirichlet boundary conditions at the 
strip [-1,1] (see e.g. [3]). 
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In pQ the solution of this equation has been represented as a series of the spheroidal 
functions which generalizes the expansion in terms of Mathieu functions for the case 
v = obtained in [4]. Equation (CQ) belongs to the class of equations whose particular 
solutions can be obtained by Latta's method [5] but the resulting equations contain 
unknown constants. For v = these constants have been explicitly calculated by Myers 
[6J through a solution of the Painleve equation of third kind (a concise summary of 
Myers' result is given in Appendix B of [7]). 

The purpose of this note is to show that methods used for v = can be generalized 
for non-zero v (|z/| < 1/2) and special solutions of ([1]) can be obtained from ordinary 
differential equations (ODE) using the Painleve III equation more general than the one 
discussed in [Bj. As a by-product we obtain an explicit formula for a solution of the 
Painleve III equation as a ratio of two infinite series of spheroidal functions with known 
coefficients. 

The plan of the paper is the following. In Section [2] we use Latta's method to 
obtain a system of ODE in t variable for two special solutions of (OQ) corresponding to 
the right-hand side of ([1]) equals coshfe and smlife. These equations contain certain 
constants which are determined in Section [3] by deriving the second system of ODE 
with respect to variable 9. The condition of the compatibility of these two systems of 
equations leads to the Painleve III equation for the constants. The limiting behaviour of 
necessary quantities for small and large 9 are derived in Section HI Though the discussed 
methods work only for the above mentioned particular choices of the right-hand side 
of ([I]), in Section [5] it is demonstrated that the knowledge of these special solutions 
permits to investigate more general problems. In Section [6] we show that kernel ([2]) is a 
positive-definite function which leads to the boundedness of the discussed Painleve III 
solution at positive arguments. For completeness, in Appendix A the series expansion of 
solutions of integral equation ([1]) is present and in Appendix B properties of spheroidal 
functions are briefly described. 



2. First system of equations 

The method of Latta [5] can be applied for integral equations of the form (CQ) when its 
kernel, K(w), obeys a differential equation whose coefficients depends linearly on w. 
The kernel K{w) determined in (T5]) obeys the equation 

wK"(w) + (1 - 2v)K\w) - w9 2 K(w) = (6) 

so the application of Latta's method is straightforward. 

Let us define the integral operator corresponding to (CQ) 

(Tg)(x) = J K(\x-t\)g(t)dt (7) 

and its two particular solutions, g c (t) and g s (t), with the right-hand side equals cosh(6 l x) 
and sinh(#x) 

(Tg c )(x) = coshes) , (Tg s )(x) = sinh(6'x) . (8) 
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It is plain that parities of g c , s {t) are fixed: g c {—t) = g c {t) and g s {—t) = —g s (t). 

An important ingredient of Latta's method is the uniqueness of the solution: when 
f(x) = the only solution of (0Q) is g(t) = 0. For kernels like (j2J) it is known that if f(x) 
belongs to L 2 (— 1, 1) then the solution is unique. A simple way to prove it is to use the 



expansion into a complete set of functions as in Appendix A 



Changing in variable w to x — t, multiplying the resulting expression by g(t), 
and integrating over t G [—1,1] gives the following identity (the first Latta equation) 

(rtg)"(x) - 9 2 {Ttg){x) = x[(Tg)"(x) - 9 2 {T g)(x)] + (1 - 2u){T g)\x) = (9) 

which permits to calculate (Ttg) when function f(x) = (T g)(x) is known and obeys the 
equation /" - 9 2 f = 0. 

The second Latta equation is obtained by remarking that the kernel in ([1]) depends 
only on the difference of arguments. Therefore, it is easy to check that for arbitrary 
function y(t) such that f/(±l) = 

(r y y(x) = (r y >)(x) . (io) 

To apply this equation it is necessary first to determine the behaviour of function g(t) 
near ends of the strip, t — ±1. Collecting the most singular terms as it is done in 
Section @] or using the representation ( 1A.5I) gives 

u 

(11) 



yC ' S t->l (1 - t)"+V2 ' 

Therefore when v < 1/2 functions 

y c>s {t) = {l-t 2 )g C)S {t) (12) 

are zero at the both ends of the strip, t/ CiS (±1) = and (flUj) is valid for them. 

Equations (El) and (fTUl) permit to derive ODE for special solutions (E]) as follows. 
By definition 

(T g c )(x) = cosh^x . (13) 
From ([9]) it follows that y = (T tg c ) obeys the equation 

y" -9 2 y = (l-2v)6smh6x (14) 

whose odd solution is 

1 - 2u 

y = (r tg c ) (x) = — - — x cosh 9x + A sinh 9x (15) 

where A is a constant independent on x. 

In exactly the same manner the knowledge of T tg c permits to calculate T t 2 g c 

(r t 2 g c ) (x) = — — — ^ — ^-x 2 cosh 9x + x sinh 9x (- — —A - - — + C cosh 9x 

8 \ 2 89 J 

(16) 
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where C is another constant. 

Similarly, with certain constants B and D 

(r g s )(x) = sinh Ox, (Ttg s )(x) — — - — x sinh Ox + B cosh Ox , (17) 

and 

(r t 2 g s )(x) = — — — ^ — —x 2 sinh Ox + x cosh Ox (- — —B - - — + D sinh Ox . 

8 \ 2 80 / 

(18) 

Using (ITUi) one gets 

(r[(i-t 2 )^')(x) = (r[(i-t 2 )^)'(x). (19) 

From the above relations it follows that 

(T[(l-f)g c ]')(x) = - (1 ~ 2z/)(3 ~ 2z/) [Ox 2 sinh Ox + 2x cosh to] 

8 

- f 1 - 2V A - 1 - ^ ) [Ox cosh + sinh Ox] + (1 - C)0 sinh 0z . (20) 
\ 2 SO / 

Expressing the right-hand side through (Tt 2 g s ), (Ttg c ) and (Tg s ) one obtains 
(r[(l-t 2 )^ c ]')(x) = -0(Tt 2 g s )(x)+ (O(B-A)-^^) (Ttg c )(x) + p x (Tg s )(x) (21) 

where pi is a certain combination of constants A, B, and C. 

As the solution of ([I]) is unique, the last equation implies that g c and g s have to 
obey the following ODE 

(1 - t 2 )g' c (t) = (i/+i- p)^ c (t) + (pi - ^ 2 )^(t) (22) 

with p = (A - B)0 + 2. 

Repeating these calculations for g s gives another equation 

(1 - t 2 )g' s (t) = {y+\ + p)tg.(t) + (p 2 - 0t 2 )g c (t) . (23) 

Imposing the condition that near the end t — 1 functions g c {t) and g s (t) have the 
prescribed singularities (cf. ( fill ) 

^ c(t) (l-t)- +1/2 ' ^ ^ (1 - 1)^ 1/2 ( ^ 
where 77 is determined from the limit 



v =tl = l im 9M (25) 



one fixes constants pi and p 2 . 
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Finally, we conclude that for \v\ < 1/2 functions g c (t) and g s (t) obey the system of 



ODE 



d_ 

dt 



9c 
9s 



M 



9c 
9^ 



(26) 



with the following 2x2 matrix M 

1 



M 



o e 
e o 



i-t 2 



t(v + l/2-p) 
V (u + l/2-p) 



~(v + l/2 + p) 
V 

t(u+ 1/2 + p) 



(27) 



When v — these equations coincide with the ones in [5]. 

For the further use we rewrite the system (T26l) in the form 



d 




t I 


dt 


9c 


r/ J 


d 






dt 


9s 


- tr] g c 



et 

— g c + 

T] 



V 



l+p 



+ 0)g 



fan 9s + (n{v - - - p) + s ) g c . 



(2f 



The terms in the square brackets are zero when t = ±1 and this system can be obtained 
directly by using (JTOjl . 

3. Second system of equations 

To find constants p and rj we use scaling arguments similar but different to the ones 
discussed in [B]. 

Functions g c and g s in (JSJ) depend on t and 0, g c>s = g CtS (t, 9). We are interesting in 
finding equations governing the evolution of these function with changing 9. 
Consider instead of system (jHJ) a more general system of equations 

TG C = / \x-t\ v K v {k\x-t\)G c {t)dt = cosh % - ±(Li + L 2 )) 

rG s = / 2 |x - *r*C(&k - t|)G s (*)dt = sinh \k(x - \{L X + L 2 ))l .(29) 
Changing the variable 

t=\(L x + L % ) + \{L 2 -L x )y (30) 

together with the similar change of x it is straighforward to check that functions G c (y) 
and G s (y) are expressed through g c (t) and defined in (jSj) as follows 



where 



L2 — Li 

L 2 + L x 



L\ L2 — i^i 



9c,s( Z iQ) 

9 = \k(L 2 — L x ) . 



Consider two functions 

My) = Gs(y) + vG c (y) , Mv) = G s {y) - vG c (y) 



(31) 
(32) 

(33) 



Painleve solution of an integral equation 

where r\ is the same as in (1251) . By construction \1/ 1 (L 1 ) = and ^/ 2 (L 2 ) = 0. 
From (EHD it follows that 



r^ lj2 = TG S ± r,YG c = sinh (k{x - \{L X + L 2 )) ± V cosh (A;(x - \{L X + L 2 )) . (34) 
As ^i(Li) = 0, the differentiation over Li of these equations gives 

d irv,) = r( ^-%) , i=i,2. (35) 



Using the uniqueness of the solutions, after simple algebra we obtain 

Performing the calculations one finds 

where i = 1 for G c and i = 2 for G s . 

Combining all terms together one obtains that fl36l) take the the form 

dg c (dg s dg, 



Using ( I2"gj) we conclude that functions g CtS {t, 0) obey the following system of equations 

(39) 

where matrix iV has the form 

1/2 + p 



TV 



, 1/2 ~P 



(40) 



The condition of compatibility of systems of equations (12~T|) and ( )39|) gives the equations 
of zero curvature for matrices N and M 

—N - —M = MN - NM . (41) 

dt dd y ' 

Direct calculations prove that this equation will be valid provided i] = rj(9) and p = p{9) 
fulfil the equations 

P^^-V'-V 2 ) , (42) 
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and 

, v + 1/2 + p , . , . 

P = " -{"+ 1/2 - P)V • 43 

V 

Substituting here the previous equation one finds that r\ has to obey equation 

dfl 2 ' W d0 6 1 r) y ' 

This equation is a particular form of Painleve III equation which has been discussed in 



From the results of Appendix A and Appendix B it follows that functions g c (cos7) 



and ^ s (cos7) can be written as the following series of the spheroidal functions 

^ oo 

£ C (C0S7) = ( sin ^y u+1 Yl ^2mA > 2m (0)y 2 m(7) , (45) 



m=0 

oo 



# s (cos7) = - — 1 )2u+1 Yl ^+1^+1(0)^+1(7) , (46) 
l sm 7 J m=0 

where F m (cos7) are defined in ( IB.lj) . X m (£) in ( IB. 181) . and /i n in ( 1A.24I) . 



Function i] in (|25|) which is a solution of the Painleve III equation is given as the 
ratio of the above series 

/X2m+1^2m+l(0)l2m+l(0) _ ^ 
Em=0 / i 2mA A 2m(0)^2m(0) 

4. Limiting behaviour 

In order to use the equations derived in the previous Sections it is necessary to know 
the behaviour of g c {t) and g s {t) defined in ([8]) in the limit of small 9 and (or) large 9. 
From §5§ it follows that at small 9 

K(w) — ► 9-"^P- + w 2 »9»~^p- . (48) 

Therefore, to find the necessary solutions of (JS]) when 9 — > it is necessary first to solve 
the equation 

\x-t\ 2 »g(t)dt = f(x) (49) 



i 



when function f(x) equals 1 or x. 

Though the general solution of such equation is known (see e.g. [12]), solutions 
with f(x) = 1 and f(x) = x can easily be obtained by the direct application of Latta's 
method [5]. After a simple algebra one finds that the solutions of equations 



i pi 



x - t\ 2v g (t)dt = 1, / \x-t\ 2u g l (t)dt = x (50) 



are 

/ \ COs(7TZ/) COs(lTu)t 
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From 0481) we obtain that 



g c (t) — > fJ, c g {t) , g 3 (t) — > l*,9i(t) (52) 
where constants // c and /x s are determined from the conditions 

P V 2 1+ ^ v^ 2 1-, T(l-z/) / ' P 2 1+ ^ v ' 

From these relations it follows that 

n(f)) __. = .1-2, cos(7n/)r(l/2 - _ 

/l ; *->o 2z//i c " ^V-^vTil - v)T(-v) 2v' K ' 

Using standard formulae for the T-function 

Y{2v) = 2 ^—Y{v)Y{v + 1/2) , T(l/2 + v)T{\/2 - u) = , (55) 

y 7T COS(7TZ/J 



one gets that the limiting behaviour of rj(9) at small 9 is 

r)(9) — > B(20) 1-2 " - — = { B ^' 1_2 "' °t" <1/2 P6) 

/v ; 0^0 v ; 2i/ 1 -0/(2v), -1/2<i/<0 V ; 

where 

p - 2-3(1-2") r (g) /c 7 \ 

^"^ r 2 (i - i/)r(2i/) ' ( j 

To find the behaviour of solutions when 9 — > 00 it is convenient to consider instead of 
functions g CjS (jHJ) a solution <7_(f) corresponding to the integral equation 

-I 

\u - v\ u K u {9\u - v\)g4v)dv = e~ eu . (58) 



-1 

To find asymptotic behaviour of g_ for large 9 we first consider the equation 

\u-v\ u K u {\u-v\)gv{v)<\v = e- u . (59) 

Its solution can be calculated either by the Wiener-Hopf method [T3] using known 
Fourier transform of the kernel 

(60) 



or by expansion formulae analogous to the ones described in Appendix A but for half-line 
integration [T5] 

nl I \x 



- t\ v K v (\\x - t\)t- v - 1 / 2 e- t / 2 L- v - 1 / 2 (t)dt 
V^r(z/ + l/2)r(n + 1/2 - u)e- x ^L- v -^ 2 {x) (61) 
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where L*(x) with n = 0, 1, . . . are the Laguerre polynomials. 
One gets 

n 

g (v) = C v v- v - 1 ' 2 b~ v , C u = \, cos ttz/ . (62) 

Using this solution it is easy to show that the dominant approximation for the solution 
of (|58p when 9 — > oo has the form 

gW(t) = C U V6(1 + t)- v - l/2 e~ te . (63) 

Straightforward transformations give 

\x - t\ u K u (6\x - t\)g {0 \t)dt = - F 2 {x) (64) 

J -i 
where 

/oo poo 
\x-t\ u K p (6\x-t\)g {Q \t)dt , F 2 (x) = J \x-t\ u K u (9\x-t\)g W (t)dt . (65) 

Changing the variable in the first integral as t = —l + v/9 and using (1591 leads to 

F 1 (x)=e" te . (66) 

In the integral for F 2 (x) it is convenient to put t = 1 + s/6. It this way one obtains 

f°° / s\ -v-l/2 

F 2 (x) = e-'-WC J (2 + -J (s + 6(l-x)yK u (s + 6(l-x))e- 9 - s ds . (67) 

When 9 is large one can drop the term s/6 in the first bracket and use the asymptotic 
form of the K u function fll]) 

POO 

F 2 (x) — ► 0-"- 1 /22-"- 1 V ^ e - JW+te a / (a + 0(1 - x)r -1/ V 2s ds . (68) 

This formula means that close to right-hand side, x = 1, the additional contribution 
related with the integration in finite limits is 

POO 

F 2 (x) « e-v- l l 2 2-»- l ^e- 2e+dx C u / s^/V^ds = e ex 6 (69) 

where 

5 = ^-9~ v - 1/2 2~ 2u cos(nis)e- 2e T(is + 1/2) . (70) 

27T 

Using fl58|) it is obvious that to cancel contribution -F 2 (x) it is necessary to modify the 
approximation ( |63i) as follows 

— ► C^[(l + t)""- 1 /^-* + 5(1 - ty u - 1/2 e te ] . (71) 

0— >co L 

This function consists of two terms. The first dominates close to t — — 1 and the second 
near t — 1. According to (|2"S"1) function 77 (0) is determined from the limit 

V (9) = ^ 9 - i r t \~ 9 - ( f\ • (72) 



Using (j7T|) we find that 



v f ) 1 _ ^!^Zr> + 1 /2)2- 2 ^^- 1 / 2 e- 2e . (73) 

6»-i>oo 7T 
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5. Embedding formulae 

In the previous Sections the integral equation ([1]) has been transformed into two systems 
of ODE fl26l) . fl27|) and fl39l) . f|40|) but only for very special right-hand sides of this 
equation, f(x) = cosh(#x) and f(x) = sinh(6'x) (cf. (j8])). Nevertheless, the knowledge 
of these two special solutions permits to find solutions of (CD) for more general cases [7j, 

Let us consider the equation of the form 

(Tg)(x) = e~ ezx (74) 

with the same kernel as in ([2]). To find its solution, g(t), it is convenient to define 
symmetric and antisymmetric combinations of g C}S 

g±{t) = g c {t)±g s {t) . (75) 

We look for solution of (I741) in the form 

*(t) = g{t) + a + g+(t)+a-g-(t) (76) 

where a + and a_ are two constants independent on t which have to be determined from 
two conditions 

Km = . (77) 
t-»±i 

The explicit form of these constants will be presented later (see (J89j)). 
Assuming that ( 1771) is fulfilled, one concludes as in Section [2] that 

(r*(aO)' = (r*')(aO ■ (78) 

The right-hand side of this equation is calculated from the definition of \l/ and 

(r*(a?))' = -6ze- ezx + 6a + e 9x - 9a„e~ dx = {Y[-z9g + a+6g+ - a_6g_]){x) . (79) 

Due to the uniqueness of the solution, ty(t) obeys the equation 

= -z0(*(t) - a +g+ {t) - a_p_(t)) + a+0g + (t) + a_6g_(t) (80) 

which is equivalent to 

rV(t) + zV(t) = (z+ l)a+g+(t) + (z- l)a-g-(t) (81) 

whose solution is straightforward if constants a± are known. 

To find them one can proceed as follows [Sj. Let us calculate the following quantities 

G7±r*) = J J g±{x)K(x-t)*(t)dtdx (82) 
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by two different methods, the first by applying the operator V on g± and the second by 
applying it on In such a manner one gets 



(83) 



J ^(t)e ±et dt = j g ± (t) (e~ zet + a+e et + a^e^dt . 
Introduce the Laplace transform of functions g CyS 

G c {p) = J g c (ty et dt , G s {p) = j g s (ty et dt . (84) 

Due to the symmetry properties of g CjS , G c (—p) = G c (p) and G s (— p) = —G s (p). 
Using (|84|) we obtain 

J V(t)e et dt = G(-z) +a+G(l) +o_G?(-l) , (85) 
®{t)e- et dt = G{z) + a+G(-l) + a_G(l) (86) 



where G(p) = G c (p) + G s (p). 

On the other hand, multiplying ( IHTi) by e ±6t , integrating the results from —1 to 1, 
and taking into account that ^(±1) = one finds 

( z ~ l )j e et m(t)dt = (z + l)a+G(l) + (z - l)a-G(-l) , (87) 

( Z + 1 )J e- et ^(t)dt = (z + l)a+G(-l) + (z- l)a_G(l) . (88) 
These and previous equations permits to calculate constants a± 

(l T z)G( T z) 

a± = — 2G(ir~ • (89) 



Combining the above formulae one finds that (181j) takes the form 

e-^'(t) + z9(t) = [G c (z)g s (t) - G s (z)g c (t)} (90) 

whose solution obeying $(±1) = is 

9{t) = y ^^) j\(y)e zey dy - G s (z) J^g s (yy ey dy] . (91) 

The above formulae relate the general plane-wave solution (|74p to two special cases 
corresponding to z = ±1. Such formulae are called embedding formulae in the theory 
of diffraction and can be derived in more general settings (see e.g. [14J and references 
therein). 
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6. Positivity relations 

Equation ( 19TT) represents the solution of ( 1T4T) . As for usual diffraction problems it is 
physically clear (and can be proved e.g. by using the series expansions as in Appendix A ) 
that such solution exist for all 9 > and any z but ( |9TI) contains G(l) in the numerator 
and in order to the solution remains finite it is necessary that this quantity is always 
non-zero. 

By the construction G(l) is the Laplace transform of the sum g c + g s 



G(l) 



(g c (t) + g s (t))e n dt 



-i 



-i 



g c (t) cosh(9t)dt + g s (t) sinh(0t)dt . (92) 



By definition, g c (t) and g s (t) are solutions of (jSJ) thus 

cosh(6 l x) = (Tg c )(x), smh(9x) = (Tg s )(x) 

and 



(93) 



G(l) = jdtj g c (t)K(\x-t\)g c (x)dx + J dt J g s (t)K(\x - t\)g s (x)dx (94) 

with the kernel K(w) defined in (J2]). An important property of such kernel is that its 
Fourier transform, K(p), presented ( 1601) is strictly positive for all p and v > —1/2. It 
means that function K(w) is a positive definite function and the double integral 



-i 



dt / f( t )K{\x-t\)f(x)dx = — 



-i 



2tt 



1 2 A 

f(t)e ipt dt K{p)dp (95) 

i 



is positive for any function f(x). Therefore for all 9 > 

> 



(96) 



and the solution fl9Tl) is finite for all z. 

The same positivity condition permits also to establish the positivity and finiteness 
of r)(9). Formally, this function is defined as the limit when t — > 1 of the ratio of two 
solutions (f25|) and it is not evident that it remains finite for all 9 > 0. 

Let us define two functions 



F e (p) = f tg c (ty et dt , F s (p) = J tg s (t)e pdt dt 



(97) 



which differ from f)84p by factor t in the integrand. 

Multiplying (|28p by e z6t and integrating the both parts over the interval [—1,1] one 

finds 



v-l/2-p 1 

9 T] 

v-l/2 + p , 1 
5 + V - 



G c {z) + -G s {z) = zF c (z) + F s {z) 
V 

G 3 (z) + z V G c (z) = F c (z) + zF s (z) 



(98) 
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From (J3HD and flU} it follows that 



^G c (z) - l!*±P-G e {z) = F s {z) + zF c (z) 



Combining these relations one concludes that 



d_ 

do 



GJz) - 



v 1 



GJz) 



rj 



-G s (z) 



d_ 

m 



GJz) 



+ T] 



(99) 



G s (z) = zr]G c (z) . (100) 



Put in the last expressions z = 1. Then G(l) = G c (l) + G s (l) obeys the equation 



v li 
9 77 



G(l) 



G(l) = C\v) 9 V exp / ( V (9') + r]-\e'))d6' 



whose solution is 



where C'(v) is a constant. 

Using d23}, d3S}, and (HU with t = 1 one concludes that 



(101) 



(102) 



k c ~\~ k c fj. 



d_, _ 1/2 + p 
Solving it and using (|42p we finds 



k 2 c =C"(u)-exp / (77(0') +V~\0')W 
V J 

with another constant C"(y). 

Comparison of this expression with f ll02p proves that 

i-f / 



(103) 



(104) 



77 fc^ = C(i/) ^-"G(l) . 



(105) 



To find the constant of proportionality, C(u), we compare the both sides of this equation 
when 0—t-O. 

From Section H] it follows that for \v\ < 1/2 and 9 — > the dominant contribution 
to G(l) is due to G C {1) and it is straightforward to check that 



C{u) 



COS 11V 



2^+l 7r 3/2 r ( 1 /2_ Z y) ' 



(106) 



For \v\ < 1/2, C(v) > and due to the positivity of G(l) it follows from ( 11051) that 
77 > for all positive > 0. As the only possible moving singularities of 77 are poles this 
reasoning proves that 77 is bounded for all 9 > (the behaviour of 77 at large 9 is fixed 
by (J73D). 
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The list of solutions of the Painleve III equation (144 j) which remain bounded as 
9 — > oo along the real axis has been presented in [10j . Such solutions form one parameter 
family parametrized by parameter A which determined the large 9 behaviour of solutions 

77(0) — > 1 - AI> + l/2)2- 2u 9- u - 1/2 e~ 2e . (107) 

Bounded solutions are characterized by the following behaviour at small 9 

7](9) — ► {29) a B + (26) B x + (29) 1+2a B 2 + {29) 2 - a B 3 + 0{9 2+a ) (108) 

6*^0 

where a (restricted by inequalities —1 < Re a < 1) is related with A as follows 

2 

a = — arcsin(7rA) . (109) 



7T 



and coefficients Bj are 



Bx = — ~ — To j B 2 = B 2 - — U -—, -, B 3 = j- -{Av 2 - (1 - of) (110) 

(1-d) 2 ' (1 + ^) 2 ^^(l-a) 41 v v ; 

with B = B(o~, v) being a function of o and v 

Bin u) - o-3 ^ 2 (( 1 -")/ 2 ) r ((l + a)/2 + ,) 

5(a ' z/) " 2 P((i + a)/2)r((i- ( x)/2 + ,) • 

When < v < 1/2 the solution of the integral equation ([I]) discussed above corresponds 
to a particular case of general theory in [10] with o = 1 — 2v (cf. ( )57|) . ( II 111) , and 
( 175]) ). It seems that the case when — l/2<z/<0is missed from the analysis of [TD] as 
1 < a = 1 — 2v < 2. For such choice of a the term i?3 in ( IllOp is identically zero and 
the asymptotics (I5^|) represents the two leading terms when 9 — > 0. Though in this case 
the dominant behaviour at small 9 is linear in 9, rj — > —9/(2i>), its asymptotics at large 
9 is given by the same formula ( 11071) but with value of A obtained by the inversion of 
(HDD (see (E3D) 

A = ^^. (112) 

7T 

7. Conclusion 

We demonstrate that special solutions of the integral equation ([T]) can be obtained 
from ordinary differential equations (|26|) . (T2T|) or (|39|) . (T4T)]) with entering constants 
calculated using the Painleve III equation (|4"4"|) with known asymptotics at small and 
large arguments. Our finding generalizes the results obtained in [6] for the special case 
v = 0. An interesting consequence of our investigation is the expression (T4T|) for a 
solution of the Painleve III equation with asymptotics given by (1561) and ( J73|) as the 
ratio of two infinite series of spheroidal functions with known coefficients. The positive- 
definiteness of the kernel forces this Painleve III solution to be bounded for all positive 
arguments thus giving another proof of connection formulae for the Painleve III equation 
[TO] in a special case. 
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Appendix A. Series solution of the integral equation 

Let ji m and Y m (j) be generalized eigenvalues and eigenf unctions of operator (JSJ) 

fi n / (sin 7 )- 2 ^(| cos/3 - cos 7 |)V m ( 7 )d 7 = Y m (f3) . (A.l) 
Jo 

As the symmetrized version of the kernel has the form 

K(P,i) = (sin P) - U K(\ cos f3 - cos 7 |) (sin 7) ~ u (A.2) 
eigenf unctions 1^(7) form an orthogonal system of functions 

(sin 7 )-^F w (7)F n (7)d7 = N m 5 mn (A.3) 



where N m is the normalization constant. 

If a function f(x) is expanded into a series of functions Y m {^) 

/(COS 7 ) = a mY m (l) (A.4) 

m 

then the formal solution of the integral equation (CQ) is 

00 

5,(0057) = (sin7) _2i/_1 ^/i m a m r m (7) . (A. 5) 

m=l 



To find explicitly eigenvalues and eigenf unctions of ( lA.lj) it was noted in [T| that the 



function K(\Jx 2 + y 2 ) plays the role of the Green function of the equation 

d 2 d 2 2v d 



dx 2 dy 2 y dy )'(■/• (A..)) 



Therefore the function 

• 1 



¥* ei \x,y) = g(t)K(^(t-x) 2 + y 2 )dt (A.7) 



is a uni-valued solution of ( 1A.6j) in all points of the (x, y)-plane except the strip [—1, 1] 



which exponentially decays at large distances 



^ (rcf) (rcos0, rsin0) — ► J^-/" 1/2 e^ C g(t)e etcos *dt . (A. 

r->oo V U7T J_ l 
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Consider a Dirichlet-type problem of finding the solution of ( )A.6j) in the form 

V(x,y) = ¥ inc) (x,y) + ¥ rei) (x,y) (A.9) 

where ^^ nc \x, y) also obeys (IA.6I) and the total field 0) is zero at the interval 
[-1,1] 

*(ar,0) = 0, -1<x<1. (A.10) 

It means that unknown function g(t) obeys ([I]) with f(x) = — ty( inc \x, 0). 

Equation ( 1A.6I) permits the separation of variables in the elliptic coordinates 



It is plain that 



X = cosh £ cos 7 , y = sinh£sin7. (A. 11) 



d 2 d 2 2 ( d 2 d 2 , , ,^ 



dx 2 dy 2 cosh2£ — cos27 V<9£ 2 d^f 



and 

I = cosh 2£ - cos 2 7 ( C ° Sh ^ Sin 7 | + Sinh * C ° S 7 Jt~) ■ (A ' 13) 
If V(x,y) = X(0^(7) then y( 7 ) and obey 

- 2z/cot 7 l^ - (a m - ^ 2 cos2 7 )F m = (A. 14) 

and 

X'^ - 2z/coth£X' 1 + (a m - i# 2 cosh20X m = . (A. 15) 

The separation constant, a m , is chosen to ensure the symmetry properties of solutions 

Y m (-j) = F m ( 7 ) , y m ( 7 ) = F m ( 7 + tt) . (A.16) 

The solution X m (£) is fixed by its normalization at infinity 

*m(0 — ► Ce^^ 1 / 2 ) exp (-^ (A.17) 

where C is a constant. 

The above equations are particular cases of spheroidal equation [9J. ^(7) is even 
periodic angular solution of (1A.14|) and X m (£) is the radial solution of the modified 
equation ( 1A.15I) decreasing at the infinity. It is known that Y m are orthogonal functions 
on interval [0, 7r] as in flA.31) and they form a complete set in the space of even functions 
on this interval. 

Therefore, the reflected uni- valued field in the (x, y)-plane cut along the interval 
[—1,1], decaying at infinity (cf. (1A.17j) ). and obeying (lA.lOj) with f(x) as in (1A.4I) can 
be represented as a formal series 



^ (ref) (e,7) = E a -f^||^)- (A.18) 
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As for the usual Dirichlet problem the value of g(t) in OA .70 is related with the normal 
derivative of \|/( ref ) at the strip [—1, 1]. At small w 

KJw) — > uT"^ + w v ^\P~ . (A. 19) 

As w = y/(x — t) 2 + y 2 and for v < 1/2 

y^o+ (x 2 + y 2 ) 1 ~ v r(l - u) V ; V ' 

one gets that [1] 

, , 2"6- u ,. o„ 9 



lim?/~ 2l/ — ^ (ref) (x,y) . (A.21) 



The behaviour of the decaying solution X m (£) at small £ follows from (IA.15I) 



X m {i) — ► A m + ,EU i+ ^ . (A.22) 
For x in the strip [—1, 1] and small y, y = £ sin 7. Therefore 

«<"^ = - ^^ (sin7) " 2 ""'£^ (7) ' (A ' 23) 

V \ / / m 

It means that eigenfunctions of operator (1A.1I) equal F m and the corresponding 
eigenvalues are 

A*™ = ~ /9 YT^ A ' 24 

0rr(l/2 - v)A m 

where A m and B m are determined from the dominant behaviour at small £ of the 
decaying solution of the radial equation X m (£) (|A.22j) . 



Appendix B. Spheroidal functions 

A convenient method to find periodic functions Y m (j) is to expand them into a series 
of Gegenbauer's polynomials [H] 

^(7)= X) kCrCcoBT) (B.l) 

n=m mod 2 

where the summation is performed over all non-negative integers of the same parity as 
m. 

Here Cfi{x) are polynomial solutions of the equation 

(1 - x 2 )y" - (2/i + + n(n + 2fj)y = (B.2) 
and C^(cos7) are solutions of 

y" + 2ficot'yy' + n{n + 2fi)y = . (B.3) 



Painleve solution of an integral equation 
We impose the standard normalization [IT] 

C£(i) 



n!r(2/x) 

so the orthogonality relation for Gegenbauer's polynomials is 



18 



(B.4) 



Notice that [EE] 



^(i _ x y wc*(x)cMx)dx - ^ n![r(/x)]2(ri + ^ 



/iC^(cos7) — > -cos(n7) 



(B.5) 



(B.6) 



so for \i — > series (IB.ip can be transformed to series used for even Mathieu functions. 
Using the recurrence relation for Gegenbauer polynomials [11 J 

(n + l)C£ +1 (x) = 2(n + n)xC£(x) - (n + 2p - l)C»_ x {x) (B.7) 

gives 

cos 2 7 ^(cos 7 ) = A n C-"(cos 7) + B n _ 2 C^ 2 (cos 7) + D n+2 (%(cos 7) (B.8) 
where 

u(l + v) „ (n-2i/ + l)(n-2i/) (n - l)n 



A 



in — v) 



, B n 



2(n-u + 2)(n-u+l) 



2{n-v- 2)(n - u - 1) 
(B.9) 

Substituting the formal series (IB.lj) into flA.14j) leads to a three-diagonals matrix for the 
determination of b n with fixed parity 



ln2 



r A, - n(n - 2u) 



ln2 



+ WB n b. 



n<Jn+2 T 2 ^nUn-2 ~ ^m^n ■ 



cv m b r 



(B.10) 



Eigenvalues of the above matrix determine separation constants a T 



From (1B.5|) it follows that the normalization constant for this solution is 

2 1 ~ 2 'Vr(n + 2/i) 



N„ 



n=m mod 2 



n![r(/i)] 2 (n + /i) ' 



(B.11) 



The function e 6 * 1 ' is a solution of (1A.6|) . Therefore it can be expanded in elliptic 
coordinates ( 1A.11I) as 

£y m ( 7 )x m (0 (B.12) 



) cosh £ cos 7 



where are certain solutions (to be discussed below) of the modified equation 

(TATT51) . 

Using (1A.3I) one gets 



X m (0 = T7- / W (sin 7 )-^e ecosh « co ^y m (7)d7 . 

M m JO 



(B.13) 
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Due to the symmetry properties 0A.160 

X 2m (0 = ^- I (sin 7 )- 2i/ cosh[0 cosh C cos 7 ]y 2m ( 7 )d7 (B.14) 

J^2m JQ 



'2m JO 

and 



X 2m +i(0 = TT— f (sin7)^sinh[ecosh^cos 7 ]y 2m+1 (7)d T . (B.15) 
Using Gegenbauer's integral [XT] 

n\ / e" COS7 C^(cos7)(sin 7 ) 2 M7 = 2' t 0Fr( y u + l/2)r(r2 + 2/i)i n ^J n+M (2) (B.16) 
Jo 



one obtains 

X m (0 = ^vi-^iA S & / (n ~ 2v) (gcosh(e))%-,(gcosh(0) (B.17) 

n=m mod 2 

where are modified Bessel functions of the first kind. 

These functions increase exponentially at large £ and the series ( IB. 171) represent a 
solution of f ]A.15j) of the first kind. As l n - v and (—l) n K n _ u obey the same recurrent 
relations, the exponentially decaying solutions called solutions of the third kind (denoted 
in 



Appendix A by X m (£)) take the form 



n=m mod 2 

(B.18) 
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